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Abstract 

Within the Improved Quantum Molecular Dynamics (ImQMD) Model incorporating the statis- 
tical decay Model, the reactions of 238 U+ 238 U at the energy of 7.0 AMeV have been studied. The 
charge, mass and excitation energy distributions of primary fragments are investigated within the 
ImQMD model and de-excitation processes of those primary fragments are described by the statis- 
tical decay model. The mass distribution of the final products in 23 8tj_|_238tj collisions is obtained 
and compared with the recent experimental data. 
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I. Introduction 



When beams in the actinide region with bombarding energies above the Coulomb bar- 
rier became available about twenty years ago, the strongly damped reactions in very heavy 
systems, such as in 238 U + 238 £7 were studied . These early experiments emphasized 

the investigation on the decay channels of the di-nuclear system (for production of super- 
heavy nuclei) or on particle creation in the strong electromagnetic fields. Recently, renewed 
interest in this subject has been motivated by the necessity of clarifying the dynamics of 
very heavy nuclear collisions at low excitation energies and by the search for new ways 
of producing neutron-rich superheavy nuclei. Based on coupled Langevin-type equations, 
a model for the simultaneous description of deep inelastic scattering, quasi-fission, fusion 
and regular fission was proposed in ref. j^]. Within this model the reactions of 238 U+ 238 U, 
232 Th+ 250 Cf and 238 U+ 248 Cm were investigated and a large transfer of charge and mass were 
found in those reactions as a result of an inverse quasi-fission process {4, 5]. Owing to very 
heavy nuclear system and very complicated process, a large number of degrees of freedom, 
such as the excitation and deformation of projectile and target, the neck formation, nucleon 
transfer, different types of separation of the composite system and nucleon emission will 
simultaneously play a role. Thus, one faces a difficulty for handling the problem with such 
complex mechanism and large number of degrees of freedom by the macroscopic dynamics 
model. In this case, a microscopic transport theory mode, is worthy to be osedflQ. In ref. 
[3] the formation and properties of the transiently formed composite systems in Strongly 
damped reactions of 238 U+ 238 U, 232 Th+ 250 Cf at £ cm =680-1880 MeV were studied based on 
the ImQMD model. One found that the weakly repulsive entrance channel potential and 
strong dissipation delay the re-separation time of a composite system, and a 15-20 MeV high 
Coulomb barrier at the surface of the single-particle potential well of the composite system 
makes the excited unbound protons still embedded in the potential well and to move in a 
common mono-single-particle potential for a period of time. These two effects restrains the 
quick decay of the composite system. That study results in our interest for the incident- 
energy dependence of lifetime of the composite system. We found that the longest average 
lifetime for the composite system of 238 U-)_ 238 u could reach to over ~1000 fm/c at the inci- 
dent energy region 1000 to 1300 MeV. Recent study on the incident-energy dependence of 
the lifetime of the transiently formed giant composite system 238 u_|_238tj by means of TDHF 
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calculations based on Skyrme energy density functional [8| confirmed this results. Since the 
correlation (fluctuation) effect is considered in the ImQMD model, it is able to calculate the 
mass (charge) distribution of primary fragments in the 238 U-|_ 238 U reaction, in addition to 
study the properties of the composite systems. The experiment for the reaction 238 u_|_ 238 tj 
at energies close to the Coulomb barrier was performed at GANIL and the mass distributions 
of products for the reaction at several energies are available now [9! flp|. which stimulates us 
to make further study of the decay of the composite system of 238 U+ 238 U. 

In this work, we study the mass distribution of products in 238 tj_|_ 238 tj a ^ 7 AMeV and 
then compare it with the experimental data. Considering the extremely complexity of the 
reaction process and saving computation time, we describe the reaction process by a two 
step model, i.e. a dynamical reaction process described by the ImQMD model followed by 
a statistical decay process which is described by a statistical decay model. 

The paper is organized as follows: In section II we will briefly introduce the theoretical 
models. In section III and IV we present the results of primary fragments and final products, 
respectively. Finally, we will give brief summary in section V. 

II. Theoretical Model 

Within this approach, the first step describes the formation and re-separation process of the 
transiently formed composite systems of 238 u+ 238 U by means of the ImQMD model. The 
primary fragments and fast particle emission are obtained at the end of ImQMD calculations. 
The second step devotes to describe the decay of the primary fragments by means of HIVAP 
incorporating with a three Gaussian model for describing the mass distribution of fission 
fragments. And finally the mass distribution of the products is obtained. 



A. The ImQMD model 



Detailed description of the ImQMD model and its applications in low energy heavy ion 
collisions can be found in Refs. jf], 7, ll|, [l2]. Here, we only mention that in this model the 
nuclear potential energy is an integration of the potential energy density functional which 
reads 

= 2^ + ^TT7T + ^ (Vp) 
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+^(p 2 -K s (V P r)6" + g T ^, (1) 
z Po Po 

where p, p n , p p are the nucleon, neutron, and proton density, 5 = (p n — p p )/(p n + Pp) is the 

isospin asymmetry. The parameters in above expressions are given in the table \l[6} . 
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TABLE I: the model parameters 

The Coulomb energy is also included in the Hamiltonian written as a sum of the direct 
and the exchange contribution: 

Ucoul = \j j p P (r)^^ Pp (r')rfrrfr' - e 2 \ " j pfdR. (2) 
In the collision term, isospin-dependent nucleon-nucleon scattering cross sections [3] are 



used and the Pauli blocking effect is treated more strictly l_4j, Il5l |. 

It is of crucial importance to make the initial nuclei in the real ground state because 
considerable excitation of initial nuclei will produce unreal particle emission and the residue 
with too high excitation which will completely masks the real decay process of residue. We 
check carefully not only the binding energy and the root-mean-square radius of the initial 
nuclei but also their time evolution. The average binding energy per nucleon of initial nuclei 
is required to be E g , s , ±0.1 MeV, where E g s , is the binding energy of nuclei in ground state. 
It is required that those initial nuclei with no spurious particle emission and their properties 
such as binding energy and root-mean-square radius being stable within 6000fm/c are taken 
to be as good initial nuclei, and then are applied in the simulation of reaction process. 
The deformation of the initial 238 U (e=0.24) is considered in the initial condition. In the 
simulation of reactions, the initial orientations of two deformed 238 U are randomly taken. 
Fig.l shows the time evolution of the binding energy and root-mean-square radius of the 
initial 238 U. 

At the end of the ImQMD calculations, fragments are constructed by means of the coales- 
cence model widely used in the QMD calculations. In this work only the primary fragments 
with mass number larger than 50 are considered. Then, we calculate the total energy of 
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FIG. 1: The time evolution of the binding energy and root-mean-square radius of the initial 238 U. 



each excited fragment in its rest frame and its excitation energy is obtained by subtracting 
the corresponding ground state energy from the total energy of the excited fragment. 



B. The statistical decay model 

The second step describes the decay process of primary fragments by emission of neutron, 
proton and a particle and fission. The statistical decay model (HIVAP code) 16| incorpo- 
rating a 3-Gaussian model for mass distribution of fission fragments for fissile nuclei is used 
to describe the decay process of primary fragments and mass distribution of final products. 
In HIVAP, the survival probability of an excited primary fragment is given by subsequent 
de-excitation process leading to a given final evaporation-residue nucleus in its ground state. 
Successive stages of a subsequent de-excitation processes for primary fragment with mass 
A, charge Z and excitation energy E are determined by branching ratios expressed by rel- 
ative partial decay widths for all possible decay modes, Ti(A, Z, E)/T to t(A, Z, E), where 
i=n,p,d,a, etc., and T to t(A, Z, E) is the sum of all particle decay widths Ti(A, Z, E) and the 
fission width V f(A, Z, E). All partial widths for emission of light particles and fission for 
excited nuclei are calculated by the HIVAP code. 

The excited actinide and transactinide nuclei in primary fragments and those produced in 
the de-excitation process undergo a fission. The production probability of a fission fragment 
with mass number A\ is calculated as follows: 



Where the P(Ai, A, Z, E) is the production probability of a fragment with mass number A\ 



from a fission of the excited nucleus with mass A, charge Z and excitation energy E. The 
P(A\, A, Z, E) is calculated based on an empirical three gaussian model. It reads 

3 

P(A U A, Z,E) = J2 9 U) (A U A, Z, E) (4) 

3=1 

and 

ff(3 ) (Al A Z E) — g^dl^g) exp[ ( Al - Aij) (A,Z,E)) 2 

g [At, A, Z, E) - ^ a(j) ^ z E ^ ex Pi 2{a U) {A ,Z,EW 1 () 

j=l,2,3. 

Where, the Gaussian distribution g^(Ax, A, Z, E) represents one of the components of 
the mass distribution of fission. Among them, the g^(Ai, A, Z, E) and g^(Ai, A, Z, E) 
describe the asymmetric component of the mass distribution, and g( 3 \Ai, A, Z, E) is for the 
symmetric component. Pw(A, Z, E), a^'(A, Z, E) and A^(A, Z, E) are the parameters for 
3 Gaussian distributions, which are the function of mass number A, charge Z and excitation 
energy E of fissile nucleus. The P^(A, Z, E) and A^(A, Z, E) obey the following relations 

P^(A,Z,E) = (1-PW(A,Z,E)) V (6) 

P^(A, Z, E) = (l- P^ 3 \A, Z, E))(l - n) (7) 
A^(A,Z,E) + A {2) {A,Z,E) = A (8) 
A^(A,Z,E) = ^ (9) 

Thus, only six parameters of P^(A, Z, E), i], A^(A, Z, E) and cr^(i=i,2,3) are inde- 
pendent, which need to be fixed according to available experimental data of fission mass 
distributions in actinide and transactinide nuclei. 

For fitting the parameters in the three Gaussian empirical formula we collect available 

n - 



experimental data of fission mass distributions [TTJ, [18|, [19j, |2Cj, |2jJ as many as possible. For 
the case of lack of experimental data the interpolation or extrapolation method is employed. 
For 238 U, data for mass distributions of fission fragments at different energies are available 
so we can obtain the energy dependence of mass distribution of fission fragments through 
interpolation. But for other fissile nuclei those data are relatively lack. For these nuclei we 
suppose that they have similar energy dependence behavior with those of 238 U because the 
corresponding theoretical study is also lack for these nuclei. This, of course, will introduce 
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FIG. 2: (Color online) The mass distributions of fission for 232 Th, 239 Pu, 252 Cf and 257 Fm. The 



experimental data are taken from 
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a considerable approximation. However, in the reaction considered in this work, the fission 
for excited 238 U is the most important one among all fissile nuclei and we expect the approx- 
imation introduced in the energy dependence of the mass distribution of fission fragments 
will not destroy the final results. In Fig. 2 and Fig. 3 we show some examples of calculated 
mass distributions of fission for different nuclei and for different excitation energies and 
make comparison with experimental data. The curves and dots are for calculated results 
and data, respectively. From the figures we can see that the empirical formula seems to be 
able to reproduce the available experimental data and able to be used to calculate the mass 
distributions of actinide and transactinide fragments. 

In order to choose matching time of two models properly, we investigate the decay 
process of the transiently formed composite systems of 238 u_)- 238 u at the energy of 7.0 
AMeV. FiG.4 shows the time dependence of the surviving probability of fragments with 
Z>110. One can see from the figure that at about t=500fm/c, two nuclei reach a touching 
configuration. After about lOOOfm/c the composite system begins to re-separate with a very 
large decay rate and at about 3000 fm/c almost all composite systems are separated. This 
process is described by the ImQMD model. The separated fragments continue to decay with 
a much smaller decay rate. This process is expected to be described by the statistical decay 
model. Thus, we select the matching time of two models to be 3000 fm/c. We have also 
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FIG. 3: (Color online) The mass distributions of fission for 238 U at different excitation energies. 
The experimental data are taken from [21] 
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FIG. 4: (Color online) The time dependence of the surviving probability of super-heavy fragments 
of Z>110. 

tried other choices such as t$= 4000fm/c, 5000fm/c and 6000fm/c and we find there is no 
change of final results. In the ImQMD calculations, 500 events per impact parameter are 
performed. 

III. The Distribution of Mass, Charge and Excitation Energy in Primary 

Fragments 

In order to study the final mass distribution of the reaction 238 U+ 238 U, we first study 
the distribution of primary fragments which are given at the end of ImQMD calculations. 
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FIG. 5: (a)The charge and (b)the mass distribution of the primary fragments of 238jj_|_238tj at 
7.0A MeV. 

The charge, mass and excitation energy distributions as well as the angular distribution 
of primary fragments are obtained by the ImQMD model calculations at time t=3000fm/c. 
The double differential cross section of a primary fragment with charge Z, mass A, excitation 
energy E and scattering angle 9 is given by: 



d 2 c P ri(Z, A, E, 
d9dE 



/ 2nbf(Z, A, E, 9, b)db = V 27r6 i A6/ i (Z, A, E, 6, b t ), 
Jo i=i 



(10) 



where fi(A, Z, E, 9, bi) is the probability of producing the primary fragments with charge Z, 
mass A, excitation energy E and scattering angle under impact parameter 6j. The maxi- 
mum impact parameter b max is taken to be 14fm since there are no inelastic scattering when 
b>14fm. The double differential cross section for primary fragments will be used as input in 
the second step for the calculations of final products in order to compare with the measure- 
ment. Let us first study the charge and mass distribution of primary fragments which is the 
integration of double differential cross sections. Fig. 5 (a) and (b) show the charge and mass 
distribution of primary fragments(A > 50) for 238 u_)_ 238 u at 7.0 AMeV, respectively. The 
sharp peak are located at the uranium for both subfigures. The primary fragments at the left 
side of the sharp peak stem from the re-separation of the composite system and fast fission 
products of actinite and transactinide fragments. The products at the right side of the sharp 
peak correspond to transuranium nuclei. The mass distributions of primary fragments at 
different impact parameter b are calculated in order to clarify the origin of the fragments 
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with different mass region. The results are shown in Fig. 6. Subfigure (a), (b), (c) and (d) 
are for the impact parameter 0-4, 5-7, 8-10 and 11-14 fm, respectively. In central collisions 
(see subfig.6(a)) the mass number distribution of primary fragments extends to A=320 with 
a big asymmetric bump around A=200-260, which means a large mass transfer between two 
uranium nuclei happening in central collisions. At semi-central collision (subfig.6(b)), the 
mass distribution becomes less wide with a much shorter tail at right hand side. There are 
two peaks appeared in the mass distribution with bigger one corresponding to uranium and 
the smaller one originating from ternary fission(occasionally from quaternary-fission)events 
in reaction 238 U+ 238 U. Clearly, in this case, very deep inelastic reaction becomes the most 
important reaction mechanism. For the peripheral collisions (subfig.6(c) and (d)) the mass 
distribution of primary fragments shows a symmetric peak with a very less variance. The 
reaction mechanism for peripheral reactions are of the inelastic or elastic scattering between 
two uranium nuclei. In order to understand the reaction mechanism and the mass distribu- 
tion of fragments evolving with impact parameters shown in Fig. 6 we present the average 
life-time of transiently formed composite system for 2 38tj_|_238tj a ^ 7AMeV as function of 
impact parameter in Fig. 7. From this figure, one can see that the life- time of compos- 
ite system increases as impact parameter decreases. In the central collisions, two uranium 
nuclei have longer interaction time, stronger dissipation of collective motion and thus have 
stronger mass transfer between them compared with larger impact parameter cases. There- 
fore, the transuranium primary fragments mainly come from the central and semi-central 
collisions. Now we study the distribution of excitation energies of excited fragments. Fig. 8 
and Fig. 9 show the excitation energy distributions for fragments with Z > 100 and 90 < Z 
< 94, respectively. As is mentioned above that the fragments with Z > 100 come from the 
large mass transfer reactions which only happen in the central and semi-central collisions, 
the results shown in Fig.8 are only for impact parameter b=0-4 and 5-7 fm. Fig.9 shows the 
results from deep inelastic scattering of 238 U+ 238 U. One sees from both Fig.8(a) and Fig.9(a) 
that the primary fragments produced in central collisions are mostly highly excited and for 
those fragments the survival probability should be very low but still there is a tail extending 
to low exciting energy, which may have certain but very small survival probability. Whereas 
for the semi-central collision (see Fig. 8(b) and Fig. 9(b)) the high excitation energy primary 
fragments decreases and the portion of low energy primary fragments increases, thus, it is 
expected that some of fragments with Z > 92 can be survival. In the peripheral collisions ( 
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FIG. 6: The mass distributions of the primary fragments at different impact parameter regions. 
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FIG. 7: The impact parameter dependence of the average life-time for the composite system of 
238 U+ 238tj at 7 AMeV 

see Fig. 9(c) and (d))the excitation energies of primary fragments are much lower compared 
with the central and semi-central cases and the reason for it is of understandable. 
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FIG. 8: The excitation energy distribution of the transuranic fragments with Z > 100 under 
condition of impact parameters (a) b<4fm and (b) b=5-7fm, respectively. 
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FIG. 9: The excitation energy distribution of the primary fragments around uranium (90<Z<94) 
for impact parameters (a) b<4fm, (b) b=5-7fm, (c) b=8-10fm and (d) b=ll-13fm, respectively. 

IV. Mass Distribution of Final Products 

From the ImQMD model calculation we obtain the distributions of charges, masses and ex- 
citation energies for all produced primary fragments in 2 38u_|_238tj collisions. These primary 
fragments will de-excite through emitting light particles, 7-rays and fission. The decay pro- 
cess and the final products are described by the statistical evaporation model (HIVAP code) 
incorporated the three Gaussian fission model described in Section II. Based on the model, 
the mass distribution of final products for 2 38tj_|_238tj a ^ ^ ne incident energy of 7.0 AMeV 
can be calculated. In Fig. 10, we show the calculated results of final products at 4 impact 
parameter regions of 0-4, 5-7, 8-10 and 11-14 fm. For central collisions( see Fig. 10(a)), the 
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re-separation primary fragments of 238 u_|_238tj systems carry high excitation energies, the 
most part of them undergoes symmetric fission and thus a single bump of mass yield is 
found at around mass number 120. The rest of fragments not undergoing fission will evapo- 
rate particles and their residues finally form a shoulder in the mass distribution around Pb, 
which is due to strong shell effect for those nuclei around Pb. The yields for transuranic 
fragments decrease rapidly as mass increase, which is due to the high excitation energy of 
primary fragments in central collisions as seen from Fig.8(a). Here we should mention that 
the yields of the transuranic nuclei is not so certain because the fission barrier and the 
fission width for super-heavy nuclei and the transuranic nuclei are largely uncertain. For 
semi-central collisions, i.e. in the impact parameter region of 5-7 fm (see Fig. 10(b)), the 
excitation energies carried by primary fragments are much less than those in the central 
collisions, so there appears a broad bump at mass number region of 80<A<170 which is 
the superposition of symmetric and asymmetric fission. There appears another small bump 
centered at Uranium(A « 230). The shallow valley between two bumps means that the 
yields of nuclei around Pb is still considerable. Here we notice that the yields of transuranic 
nuclei is relatively higher than those in central collisions, which is because the excitation 
energies of primary fragments are much lower than those in central collisions. For peripheral 
collisions (see Fig. 10(c) and (d)), elastic or inelastic scattering play a dominant role and the 
behavior of low energy fission of actinide nuclei is shown. The small shoulder around Pb 
seems to appear for impact parameters b=8-10fm(see Fig. 10(c)). 

In order to make comparison with experimental measurement we have to make selection 
of scattering angle to fit the angle cut in experimental data, i.e. only fragments with the 
scattering angles of 56° < 6 < 84° and 96° < 9 < 124° in the center of mass frame are 
selected[9j. In the calculations, we assume that the scattering angle of residue of the primary 
fragment which undergoes emission of light charged particles is the same as the fragment 
itself. This assumption is roughly reasonable since the mass of residue is much larger than 
that of emitted light particles. For fragments from fission, we assume that the outgoing 
angle of one fragment is randomly distributed in the rest frame of the fissioning nucleus and 
the outgoing angle of the other one is then obtained by momentum conservation. Finally, we 
obtain the mass distribution of the final products with the same scattering angle cut as that 
in the experiment. The results are shown by open triangles in Fig. 11. The experimental 
mass spectra from [o| are also indicated by solid squares, open squares, solid circles, open 
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FIG. 10: The mass distribution of the products in 238 U+ 238 U at different impact parameter 
regions. 

circles and solid triangles for incident energies of 6.09, 6.49, 6.91, 7.10 and 7.35 AMeV, 
respectively, in Fig. 11. From the figure we find that the behavior of the calculated mass 
distribution at 7.0 AMeV is generally in agreement with the data at incident energy 7.10 
AMeV except the yields at mass region from 170 to 210 to be overestimated compared 
with experimental data. The following most important features of mass distribution can be 
deduced: (1) A dominant peak around uranium is observed, which can be attributed to the 
contribution of the reactions with large impact parameters seeing from Fig. 10; (2) The steep 
decreasing yield above U with the increase of the mass number is appeared. The products 
at this mass region stem from large mass transfer in small impact parameter reactions; 
(3)A small shoulder can be seen in the distribution of the products around Pb, compared 
with the products with mass near and smaller than Uranium for which the yields decrease 
exponentially as mass decreases. The appearing of small shoulder around Pb is due to the 
very high fission barrier around Pb. The central, semi-central collisions and even reactions 
with b=8-10fm contribute to the shoulder around Pb region; (4) In the region below A^190, 
the double bump distribution is observed. These products are from the fission of actinide 
and transuranium nuclei obviously, which are superposition of symmetric and asymmetric 
fission. 

V. SUMMARY 

In this paper, we apply the microscopic transport model, namely the ImQMD model 
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FIG. 11: (Color online) The mass distribution of the products of reaction 238 U+ 238 U. The experi- 
mental data are from [{J. 

incorporating the statistical decay model and empirical three Gaussian fission model to study 
the reaction mechanism and the mass distribution of products in the reaction 238 U+ 238 U at 
the incident energy of 7.0 AMeV. The mass, charge and excitation energy distributions of 
primary fragments are calculated within the ImQMD model and the de-excitation process 
of those primary fragments is studied by using the statistical-evaporation model (HIVAP 
code). The impact parameter dependence of the mass distribution of primary fragments and 
final products are analyzed, from which the origin of products at different mass region can 
be understood. Finally, the mass distribution of final products in 2 38tj_|_238tj collisions with 
scattering angle cut is calculated in the first time and compared with recent experimental 
data. The main features of experimental mass distribution are reproduced, those are: (1)A 
dominant peak around the uranium nuclei is observed, which corresponds to elastic and 
quasi-elastic reaction products; (2) The yields of the transuranium nuclei decrease rapidly 
with increase of the mass A ; (3)A small shoulder can be seen in the mass distribution 
of the products around Pb on the background of products for which their yields decrease 
with their mass deviating from uranium exponentially. Those products are the residues 
of primary fragments surviving from multiple-particle evaporation; (4) In the region below 
A^200, the double hump mass distribution is observed, which are the fission products from 
superposition of symmetric and asymmetric fission and mainly come from the fission of nuclei 
around uranium and transuranium fragments at high and low excitation energies. The main 
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discrepancy of our calculated results with experimental data is overestimate of mass yields 
in region of 170-200 and underestimate of mass yield of transuranium nuclei, which mainly 
come from the calculation of fission mass distribution for actinide and transuranium nuclei 
and fission width at high excitation energies. Further study is still under a way. 
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